Bandgap renormalization and excitonic binding in T-shaped quantum wires 
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We calculate the electronic structure for a modulation doped and gated T-shaped quantum wire 
using density functional theory. We calculate the bandgap renormalization as a function of the den- 
sity of conduction band electrons, induced by the donor layer and/or the gate, for the translationally 
invariant wire, incorporating all growth and geometric properties of the structure completely. We 
show that most of the bandgap renormalization arises from exchange-correlation effects, but that a 
small shift also results from the difference of wave function evolution between electrons and holes. 
We calculate the binding energy of excitons in the wire, which breaks translational invariance, using 
a simpler, cylindrical model of the wire. For a single hole and a one dimensional electron gas of den- 
sity n e , screening of the exciton binding energy is shown to approximately compensate for bandgap 
renormalization, suggesting that the recombination energy remains approximately constant with n e , 
in agreement with experiment. We find that the nature of screening, as treated within our non-linear 
model, is significantly different from that of the various linear screening treatments, and the orthog- 
onality of free carrier states with the bound electron states has a profound effect on the screening 
charge. We find that the electron and hole remain bound for all densities up to ~ 3 x 10 6 cmT x and 
that, as n e increases from zero, trion and even "quadron" formation becomes allowed. 
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I. INTRODUCTION 



The gap in a semiconductor heterostructure between 
conduction and valence bands, and the interaction be- 
tween electrons in the one and holes in the other, are 
known to depend in a complicated fashion on the pres- 
ence of mobile charges |l]-|3|] . The interpretation of optical 
experiments in doped semiconductor quantum wires , 
for example, must invoke the variation of band edges due 
to many-body and geometrical effects, the "red shift" 
associated with electron-hole binding into excitons and 
finally the reduction of the exciton binding energy due 
to screening by free carriers. 

Theoretical description of bandgap renormalization 
(BGR) and exciton formation and screening are fre- 
quently addressed with a many-body formalism HQ. 
Within this framework, contribution of both electron- 
phonon and electron-electron self-energies to BGR || as 
well as the influence of dynamical screening on the ex- 
citon binding energy can be studied. The many-body 
treatment, however, has the disadvantage that, for BGR 
it commonly ignores geometrical factors, such as the 
quantum confined Stark effect whose relevance is 
structure-specific. Furthermore, in the exciton problem, 
many-body theory treats screening within the linear ap- 
proximation and, generally, influence of the bound elec- 
tron on the free electrons is not fully included. In par- 
ticular the orthogonality of the free electron states with 
the bound state, which increases in importance in lower 
dimensional systems, are typically not included JTc[| . 

A different formalism without these shortcomings, al- 



beit one which forsakes the electron-phonon interaction 
and the dynamics of screening, is that of density func- 
tional theory (DFT) Q. Within DFT the ground state 
energy of an interacting many particle system is known 
rigorously to be given by a functional of the density. 
The essential problem of DFT is that this functional 
is unknown and the components beyond kinetic and di- 
rect electrostatic energy are isolated into an "exchange- 
correlation" functional which must be approximated rad- 
ically. Typically the limit wherein the density varies 
adiabatically, known as the local density approximation 
(LDA), is assumed. Nonetheless, successful treatments 
of a wide class of systems; from atoms and molecules 
to solids and heterostructures, abound in the literature. 
It is the purpose of this paper to provide a theoretical 
description of BGR and exciton screening, applied par- 
ticularly to semiconductor quantum wires (QWs), within 
DFT. 

We are interested in quantitative comparison with op- 
tical experiments on T-shaped quantum wires (TQWs) 
in the presence of a one-dimensional gas of electrons 
(1DEG) with ID density n e , induced in the conduction 
band via a combination of modulation doping and gating 
Jl2| . By way of warning, this one-component plasma con- 
trasts with the case of many studies which focus on intrin- 
sic QWs wherein an overall charge neutral electron-hole 
plasma is generated entirely through photoexcitation. In 
our case photoexcitation is assumed to provide a small 
number of holes (which we take as a single hole) and to 
have a negligible effect on the density of conduction band 
electrons. 
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FIG. 1. Schematic of T-shaped quantum wire cross section. 

We evaluate the band structure of a TQW in two steps. 
In the first step (section II) we consider a realistic model 
of a cross-section of a TQW, shown in Fig. 1, and solve 
self-consistently for the electronic structure of the sub- 
bands, as a function of the gate voltage, assuming trans- 
lational invariance along the wire. Material-specific ef- 
fective masses and band offsets, as well as the actual 
growth profiles and doping densities of the structure are 
included. The exciton problem, which breaks translation 
invariance, is treated in the second step (section III) by 
a simplified, cylindrical model of a quantum wire which 
employs the results of the first step for its rationale and 
its parameters. Thus the philosophy of the calculation is 
that the exciton formation and its screening represent an 
additive correction to the translationally invariant band 
structure of step 1. This exciton effect can be analyzed 
within a single subband, generic wire model so long as 
the effects of lateral confinement are correctly borrowed 
from the more faithful, albeit 2D calculation in step 1. 

Among our main results are the following. For the 
translationally invariant band structure (section II): (1) 
We find a BGR as a function of n e with comparable order 
of magnitude to that derived within many-body theory. 
(2) One small (~ 20 %) component of this BGR arises 
from a variation of the lowest electron subband eigenfunc- 
tion with n e . Specifically, at higher n e the lowest elec- 
tronic wave function spreads somewhat along the over- 
growth interface (throughout QW2 in Fig. 1) whereas it 



is more strongly confined at the junction (i.e. the inter- 
section of the two wells) at lower n e . The hole wave func- 
tion, by contrast, is always localized at the junction. (3) 
The stability of the hole eigenstate at the junction is de- 
pendent on the boundary conditions used for large z (see 
Fig. 1). If the chemical potential is pinned to an assumed 
background donor density in the substrate (we use the 
term "substrate" to denote the region far from the donor 
plane, i.e. large z, see Fig. 1), then the hole bound state 
becomes a resonance which is metastable to escape into 
the substrate. On the other hand, when the substrate is 
treated as exactly electrically neutral (or pinned to ac- 
ceptors) the hole state at the junction is stable. (4) For 
the specific dimensions and compositions which we con- 
sider |jl2"f the wires are predominantly one-dimensional 
with filling of no more than two subbands and spreading 
of the density into QW2 generally less than ~ 15 % of 
the total density. 

Regarding the exciton problem (section III), we find: 
(1) the exciton binding energy weakens considerably as 
the 1DEG density is turned on, but never vanishes com- 
pletely. (2) The requirement of orthogonality between 
the bound electron eigenfunction and those of the screen- 
ing electrons leads to an anti- screening of the exciton, 
such that the density of free electrons in the vicinity of 
the hole is depressed. This effectively positive charge 
raises the energy of the hole while the energy of the elec- 
tron actually becomes more negative and its localization 
about the hole increases with density. 

Finally, the combination of the bandgap renormaliza- 
tion and the exciton screening tend to cancel and, consis- 
tent with experiments, the recombination energy remains 
relatively constant with n e . This result is plotted in fig- 
ure 6 but discussed at the end of section III. 

We principally employ effective atomic units wherein 
1 Ry* = m*e 4 /2h 2 K 2 » 5.25 meV and 1 a* B = 
h 2 n/m*e 2 sa 100 A. For comparison with experiments 
we also use cm -1 , where 1 a*^ 1 ~ 10 6 cmT 1 . 

II. TRANSLATIONALLY INVARIANT WIRE: 
BGR 

We calculate the electronic structure of the GaAs — 
AlGaAs TQWs [[l3| by solving Schrodinger and Poisson 
equations self-consistently, within a region illustrated in 
Fig. 1, for the conduction band electrons, and by in- 
cluding exchange and correlation in the local density ap- 
proximation (LDA) (ll]]. There are numerous calcula- 
tions of BGR in lower dimensional systems which em- 
ploy many-body theory to determine the self-energy cor- 
rection to the subbands from exchange and correlation 
effects 0,^5) • In DFT a bandgap renormalization arises 
as follows. First, variations in the electrostatic potential 
tend to affect both bands equally except for kinematic ef- 
fects like the quantum confined Stark effect M where the 
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electron wave function is centered on a different spatial 
location (at different potential) than that of the hole. By 
contrast, exchange and correlation effects lower the en- 
ergy of electrons in proportion to n e , the electron density, 
and raise the energy of holes in proportion to n^, the hole 
density, producing an overall shrinkage of the bandgap. 
For the one component plasma considered here, there- 
fore, the exchange-correlation potential acts only on the 
conduction band electrons, since rih ~ 0. 

One advantage of our procedure is that while many- 
body calculations in principle begin from wave functions 
and subband eigenenergies derived from some Hartree 
calculation, the variation of these properties with den- 
sity and the aforementioned kincmatical corrections to 
the bandgap are usually ignored. By contrast DFT auto- 
matically includes these. Furthermore, the separate con- 
tributions from kinematic effects as opposed to exchange 
and correlation are easily isolated (to first order) by per- 
forming a calculation without the exchange-correlation 
potential (i.e. Hartree only). 

We employ the parameterization of Ceperley and Alder 
p6| for the density-dependent, exchange-correlation po- 
tential V xc (n(x.)) . Taking y as the translationally invari- 
ant direction (see figure 1), the Schrodinger equation in 
the x — z plane (k y — 0) reads 
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-V + V c (x,z) + e(/)(x,z) + 



2m*(x, z) 
V xc {n(x,z))]ip n (x,z) = E^Tp n (x,z) 



(1) 



where Vc(x,z) is the conduction band offset |17j] which 
depends only on the local aluminum concentration r\ (we 
use the symbol r\ instead of the standard x to avoid con- 
fusion with the coordinate x), the 2D eigenfunctions and 
eigenvalues are given by ip n (x,z) and E n respectively, 
and the electrostatic potential is <p{x,z). Also in equa- 
tion [j] the level index is n, and the subscript e refers 
to electrons in the conduction band whereas for holes we 
use E?, The eigenstates are abbreviated below as en and 
hn for electrons and holes, respectively. Periodic bound- 
ary conditions are taken at x = ±a/2 where a = 5 a' B 
is the period of the TQW superlattice as fabricated [12| . 
The effective mass dependence on position arises through 
its dependence on r\ JT^ ]. In the specific device we con- 
sider, all regions have either n = (pure GaAs) or else 
rj = 0.3 for the barrier regions. We ignore interface grad- 
ing effects as well as image effects, taking the dielectric 
constant as k = 12.5 everywhere. 

For the solution of Poisson's Equ. we also assume pe- 
riodic boundary conditions at the x-borders. For the 
surface z — we simulate the surface metal gate in the 
standard fashion |L3| by Dirichlet boundary conditions 
which fix the gate potential modulo an offset of 0.8 eV 
for the Schottky barrier. Variation of the electron sheet 
density in the overgrowth well N e (or, equivalently, the 
ID density n e ) is accomplished by biasing this surface 



gate. Other experimental methods for varying n e include 
illumination or a variation of the thickness of the spacer 
between the modulation doping layer and the overgrowth 
well. While we have modeled each of these methods inde- 
pendently, in fact in each case there is an approximately 
homogeneous sheet (or sheets) of positive charge balanc- 
ing the electronic charge in the wire/ well and the dif- 
ference in the electronic structure between the various 
methods, for a given resultant n e , is negligible. Thus 
here we actually vary a surface gate voltage to change 
density, but in displaying the results it is sufficient just 
to plot variables versus n e . 

For large z we find, as noted, that the solution for 
the holes is sensitive to the asymptotic value of the po- 
tential. Generally, we assume Neumann boundary con- 
ditions, which is equivalent to assuming complete neu- 
trality of the substrate. However we discuss below the 
nature of the hole metastable state in the case where the 
chemical potential is pinned to a shallow donor level at 
large z. 

Since the modulation doping is n-type, the Fermi level 
is close to the conduction band throughout the device. 
Furthermore, since even during photoluminescence mea- 
surements, the excitation power is very low, we ignore the 
valence band contribution to the charge density. There- 
fore, the total density which enters the Poisson equation 
is given as 



p(x, z) = N+(x, z) + n(x, z) 



(2) 



where n(x, z) is the density of conduction band electrons 
only, and where the charge density of ionized donors is 
N + (x, z). The conduction band electron density is given 
in the well regions by 



n(x, z) 



j dk Y,f^-kl-E n )\^ n {x,zt 



(3) 



where / is the Fermi function, fi is the chemical poten- 
tial (assumed constant throughout the device and taken 
as the zero of energy) and ky is the kinetic energy of 
motion in the y-direction. We also employ a Thomas- 
Fermi approximation for the density of electrons in the 
barrier regions away from the wells. Finally, we assume 
temperature T — 4 K . 

We find little difference in the results whether the 
donors, which are silicon atom DX centers |fl8| , are 
placed in thermal equilibrium with the electron gas or 
are treated as a simple sheet of positive charge with fixed 
areal density. Insofar as at low temperatures the latter 
approximation is more physical, we use this assumption 
in the calculations discussed in this paper. 

Upon obtaining a self-consistent solution of 
Schrodinger and Poisson equations for conduction band 
electrons, we calculate the hole eigenfunctions with the 
potential 



Vh(x, z) = Vv(x, z) - e<f>(x, z) 



(4) 



3 



where now <f>(x, z) is the electrostatic potential which 
solves the self-consistent problem above. The valence 
band offset Vv(x, z) also depends on the aluminum con- 
centration. We assume a bandgap offset parameter Q e = 
0.6. We treat the holes in the simplest approximation, 
assuming a single, heavy-hole isotropic effective mass 
which depends only on aluminum concentration; specif- 
ically m* h {GaAs) = 0.377 mo and m^(Alo. sGclo.tAs) = 
0.403 m H|. 

To summarize therefore and to make the context of 
this calculation clear, we are interested principally in the 
BGR as a function of n e , the integral over x and z of 
n(x, z), which we imagine to be modulated experimen- 
tally with an electrostatic gate on the overgrowth sur- 
face @. Thus we solve the evolving self-consistent elec- 
tronic structure assuming the chemical potential to be far 
from the valence band (i.e. no holes). We solve for the 
hole subbands as a one particle problem after solving the 
electronic structure at given n e , using the resultant elec- 
trostatic potential. We use a simplified band structure 
for the holes, in contrast to ref. |^0| which employs the 
full, four-band Luttinger Hamiltonian but which assumes 
n P = 0. 



A. Electronic structure 

The parameters for the calculation are summarized in 
the schematic, Fig. 1. We follow reference by desig- 
nating the overgrowth quantum well as "QW2" and the 
initial growth well as "QW1." 
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FIG. 2. Subband energies vs. wire ID density for electron 
(e) and hole (h). Electron levels measured with respect to 
Fermi surface; hole level with respect to escape to substrate. 
Electron continuum begins at dashed line. Filling of El be- 
gins at n e ~ 0.5. 

The structure of the subbands of conduction band elec- 



trons as a function of n e is shown in Fig. 2. The zero 
of energy for the electrons is the Fermi energy Ep (the 
plotted hole energy is discussed further below). In the 
higher density regime a second subband, which has con- 
centration in the inter- wire QW2 region (cf. Fig. 4), is 
occupied. The region marked "continuum" corresponds 
to the beginning of a dense set of states which are asymp- 
totically free to escape into the substrate (i.e. along QW1 
which runs in +z direction). This is so because QW1 is 
wider than the QW2, and the electrostatic advantage of 
proximity to the positive dopants in QW2 is overcome 
by the additional confinement energy there. Nonethe- 
less, the spreading of the occupied states, when the total 
density is increased, occurs within QW2 (so as to com- 
pensate the positive charge of the donors and gate). 

It is noteworthy that for n e ~ 0.5 a*^ 1 (« 0.5 x 
10 6 cm -1 ), there are two subbands which are below the 
continuum in contrast to the situation for lower n e where 
only a single bound state exists. This is also consis- 
tent with the bare (n e = 0) TQW studied in Ref. |2C]] 
which had a nearly identical aspect ratio (i.e. ratio of 
well widths) and exhibited only a single state separate 
from the continuum. The estimate there of a separation 
between E® and continuum of ~ 15 meV is also consis- 
tent with our slightly lower value of ~ 13 meV. Inter- 
estingly, this separation remains approximately constant 
(Fig. 2) as n e is increased, until E\ drops down from the 
continuum and begins to fill. 
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FIG. 3. Areal density (integrated along z) vs. x for various 
total wire densities. Confinement to wire almost complete 
below n e ~ 0.6. For highest n e interwire density develops 
subsidiary maximum at ±2.5 a* B . 

In Fig. 3 we have plotted the 2d surface density as 
a function of x for various total n e . As n e approaches 
1 Ojf 1 some spreading of the density away from the T- 
junction occurs. This results both from the occupation 
of el and also from slight wave function spreading of eO 
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(cf. Fig. 4). Note, however, that, rather than spreading 
smoothly out from the wire, the density achieves a local 
maximum at an inter-wire minimum in the potential (at 
±2.5 a* B which, due to periodicity, are equivalent points), 
analogous to the electrostatic potential minima in the 
barriers of a semiconductor superlattice. Here, however, 
the effect is both electrostatic and quantum mechanical. 
Note, in figure 4, where the moduli squared of the eigen- 
functions are plotted, that the subband el has two nodes 
along the x-direction (i.e. along QW2). For different pa- 
rameters (not shown), such as for a wider QW2 or much 
higher n e , a state with a single node at the T-junction 
is lower in energy. But in this regime the attractiveness 
of the T-junction is sufficient to stabilize the even node 
state. 




FIG. 4. Wave functions (moduli squared) showing spread- 
ing of eO as n e increases from 0.12 x 10 6 cm -1 (LR) to 
0.9 x 10 6 cm -1 (LL). First excited electron state at shows 
two nodes along x (UL). Hole ground state (UR) unchanging 
and strongly confined at junction. 

The spreading of eO with increasing n e observable in 
Fig. 4 is relevant to the evolution of the gap between 
conduction and valence bands. We will see below that 
the dominant portion of the BGR arises from exchange 
and correlation effects. Nonetheless, as noted in the in- 
troduction, the difference between the evolution of elec- 
tron and hole wavefunctions with n e produces a kind of 
kinematic BGR as in the quantum confined Stark effect. 
Explicitly, the hole eigenfunction, which is shown in Fig. 



B. Hole states 



The stability of the hole state at the T-junction de- 
pends weakly on the chosen boundary conditions for Pois- 
son Eq. at large z. If we pin the chemical potential to 
an assumed shallow donor level at large z, a common 
assumption, the hole state at the T-junction is merely 
metastable, and escape, through a long, shallow bar- 
rier, to the substrate is energetically preferred. Assum- 
ing Neumann boundary conditions, or "flat bands" at 
large z, on the other hand, results in a true bound state. 
Dynamically, the difference between "pinned band" and 
flat band conditions is mostly inconsequential for holes 
photo-generated at the junction. For any reasonable 
background (i.e. un-intentional) donor density, the bar- 
rier to escape of the hole is too large and the hole will 
remain and, presumably, recombine in the wire. Further- 
more, the hole will induce an image in the electron gas 
which will cause further, electrostatic binding to the wire 
region. A similar effect of holes bound to a 2DEG has 
been discussed recently in Ref. |^l| . 

For holes generated in the substrate, far from the wire, 
however, the dynamics of the diffusion of holes into the 
wire region can be realistically expected to depend on the 
band shape as reflected in this dichotomy over boundary 
conditions. Generally, fewer background donors is favor- 
able to hole diffusion from substrate into wire. 




4 only for n e = 0.9 x 10 6 



FIG. 5. Contour plot of electrostatic potential for 
undergoes essentially n e = 0.9 x 10 6 cm' 1 . Schematic illustrates field lines, which 



no change with n e over the n e range of Fig. 2. There- 
fore the hole essentially tracks the electrostatic potential 
at the T-junction. The electrons in eO however, for in- 
creasing n ei are able to lower their energy by spreading 
along QW2. This implies a redshift in the bandgap which 
emerges even in the absence of exchange-correlation ef- 
fects (see the Hartree results in Fig. 6 below). 



help to localize hole at junction. 

Even for pinned band conditions the basic charge dis- 
tribution near the wire establishes a purely electrostatic 
barrier to hole escape, as seen in Fig. 5. Here we plot 
the 2d electrostatic potential for n e = 0.9 x 10 6 cm -1 . 
The potential contour principally derives from a plane of 
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positive charge (the donor layer) and the line of negative 
charge in the wire. This results in an electrostatic poten- 
tial hill near the junction (to which holes are attracted). 
Thus the hole has at least a long shallow barrier to escape 
even before considering the quantum effect that further 
binds the hole to the junction. 

In Fig. 2 we have plotted the energy of the first hole 
subband, E®, versus n e . In this case we have used Neu- 
mann boundary conditions for Poisson's Eq. at large z 
and the energy zero for E® is the electrostatic potential 
at large z added to the confinement energy of QW1; in 
other words it is the energy required to escape to the sub- 
strate continuum. Here we see that the state is stable, 
and not merely metastable. Further, the increase in the 
binding with n e (up to n e ~ 0.6 x 10 6 cm -1 ) arises from 
the strengthening of the electrostatic potential minimum 
(for holes) shown in Fig. 5. 



C. Bandgap 




Density n e (10 cm ) 

FIG. 6. Bandgap renormalization for full LDA (solid) and 
for pure Hartree (dashed), here in meV, as a function of den- 
sity. Exciton binding energies for wire radius a = 0.55 a* B 
(pluses) and a = 0.35 a* B (crosses) and LDA bandgap cor- 
rected for excitonic energy as squares and triangles, respec- 
tively. 



The bandgap in the heterostructure depends upon the 
assumed intrinsic bandgap in the host materials, GaAs 
and AlGaAs which, as noted, we have taken from the lit- 
erature |Tij] . The bandgap variation with n e is obtained 
from a sum of the lowest subband energies for the elec- 
tron and the hole, each measured with respect to the 
GaAs band edge at any common point, then added to 
the intrinsic bandgap. Figure 6 illustrates the variation 
of the bandgap as a function of n e . The rapid change 
for small n e for the full LDA calculation derives from 
the form of the exchange-correlation potential. The ex- 
change energy alone, for comparison, is proportional to 
the Fermi momentum kp — sfEp p2[ . The overall or- 
der of magnitude and shape of the BGR are consistent 
with the results of calculations within the GW approx- 
imation [|| . However, since the authors of reference ||] 
assume a neutral, two-component plasma, the results are 
not quantitatively comparable. 

Also in figure 6 we plot the BGR for the case where the 
exchange-correlation potential is turned off, that is, for 
pure Hartree. Here the bandgap still varies, albeit much 
less, with n e . The origin of the shift, alluded to in the 
discussion of Fig. 4, is the variation of the wavefunction 
of the electron with n e which causes a kinematic compo- 
nent to its energy while the hole eigenstate is essentially 
dormant at the junction. This effect is a slight general- 
ization of the quantum confined Stark effect, wherein the 
bandgap is affected by an electric field, which is varied, 
between the electron location and that of the hole. In 
our case the electron additionally shifts its location as 
the gate bias varies. 



We postpone a discussion of the final result of Fig. 6, 
the behavior of the BGR when correction for excitonic 
energy and screening are included, until after the next 
section. 



III. EXCITON 

The electronic structure of the exciton can be treated 
in a fashion analogous to the calculation of electron bind- 
ing to shallow impurities, which was described some time 
ago by Vinter |23| for a silicon inversion layer. In con- 
trast to the impurity problem though, the hole cannot 
be treated as a point charge, but has a wave function 
which must be computed self-consistently with those of 
the bound and free electrons. 

In this section we employ a simplified geometry, shown 
in Fig. 7, for a single mode, cylindrical quantum wire 
along direction y of radius a surrounded by a cylindrical 
metal gate at radius R whose potential U ex t is used to 
vary the electron density. The radial charge distribution 
of both electrons and holes is assumed to be given by a 

single Gaussian wave function £j (r) = \ \J~^&~ r ^ a ■ The 
eigenstates of the electrons are determined by 

(-^^2 + v*ff(v) + v xc (v))Mv) = £ p ei> P (y), (5) 

where 

V eff (y) = J d 2 v\^{r)\ 2 eUr,vl (6) 
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and where the electron effective mass is m e and <f> e (r, y) is 
the electrostatic potential (see below) . The length of the 
cylinder is L y and Eq. || is treated with periodic bound- 
ary conditions at the ends y = ±L y /2. The effective ID 
exchange correlation potential at y is determined from 
the 3D parameterized form by 



V xc (y) = 2Tr / drrrp e (r,y)V xc (p e (r,y)) 



(7) 



where p e (r,y) — \(,i(r)\ 2 p e (y) and where p e (y) is deter- 
mined by filling eigenstates of Eq. ^| up to a fixed chemi- 
cal potential p, taken as the energy zero of the problem. 
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FIG. 7. Schematic of cylindrical wire for model exciton cal- 
culation. Hole centered at midpoint of wire. Radial charge 
distribution of hole and electron taken as single Gaussian. 
Drawing not to scale: b > 10a. 

The electrostatic potential in || is determined by solv- 
ing Poisson's equation 



4-7T 

V 2 e (r,?/) = — P(r,y) 

K 



(8) 



subject to Dirichlet boundary conditions at the gate, 
r = R. and Neumann conditions at y = ±L y /2. Here 
p(r, y) is the charge density of both electrons and the 
single hole. 

An expression similar to Eq. without the exchange- 
correlation potential, and with the hole mass in place 
of m e , can be written down for the hole, where only the 
lowest eigenstate is filled. Note however that in order 
to avoid self-interaction of the hole, the electrostatic po- 
tential (and the y-dependent effective potential) in this 
case is calculated from p e (r,y), i.e. the density of only 



the conduction band electrons, and is denoted cj>h- Self- 
interaction of the electrons is assumed to be compensated 
for by the exchange-correlation potential, as is standard 
in DFT; although it is known that the cancellation in 
LDA, unlike Hartree-Fock, is not exact. Note that a 
small inaccuracy of the model is that, by using, for the 
hole, the electrostatic potential computed without the 
hole charge Ph(r) included, the image charge induced by 
the hole on the gate, with which the hole physically does 
interact, is absent. A more precise, though computation- 
ally taxing, approach would be to compute 4>h(r, y) from 
all the charge but then subtract J d 3 r / Pfe(r') | r _ r / | i tne 
potential of the hole in free space. 

The binding energy Eb of the exciton is defined as the 
difference in the total energy between the case where the 
hole is localized, on the one hand, and the case where 
the single hole charge is spread uniformly in y, a free 
hole, on the other hand. The total energy must include, 
in addition to a sum of occupied eigenvalues, a double- 
counting correction as well as the energy related to the 
gate. Furthermore, this gate energy is modified by the 
work supplied to the gate [{13 24 1, and is hence really a 
free energy. The overall expression for this free energy is 
thus 



F(U ext ) = el + Y, n P £P e 



I 



d 3 r( Pe (r)<t> e (r) + Ph (r)<t> h (r)) 



2 . 

+ J d 3 rp e (r)(e xc {p e (r)) -V xc {p e {r))) 



(9) 



where Q is the total charge induced on the gate, deter- 
mined by computing the normal derivative of the poten- 
tial at the gate surface, and where n p is the occupancy of 
state p given by a Fermi function. Also, the third term 
on the right hand side of Eq. || is a form of double count- 
ing correction for the exchange-correlation energy, where 
£xc{p) is the exchange-correlation energy per particle of 
a homogeneous electron gas of density p |fl| . 

We note that one might write the total energy as sim- 
ply the sum of occupied Kohn-Sham energies (including 
the hole) via Koopman's theorem. To lowest order, where 
only a single bound electron and the hole energies are 
changing appreciably, one could further estimate this as 
simply F(U ext ) = e\ + e\. However these are both ap- 
proximations and the correct formula, Eq. must be 
used for all but heuristic purposes. 



A. Results 

In figure 8 we show Eb computed as a function of U ex t 
for a = 0.35, 0.45 and 0.55 a* B - Also plotted is the 
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FIG. 8. Binding energy of exciton vs. gate potential for 



this shift occurs, corresponding to a single electron in 
the wire. Clearly the minimum depends somewhat sen- 
sitively on the choice of a, and the binding can become 
quite strong for very narrow wires, as is well known [p5[ . 
The small fluctuations visible in the energies result from 
the finite size L y and the discreteness of the "free" states. 

Of considerable interest is the absence of a binding to 
unbinding transition for the exciton. Even when we ex- 
tend the density range to > 2 a'g 1 (i.e. > 2 x 10 6 cmT 1 ) 
and even for the widest charge distribution, a = 0.55 a* B , 
the ground state remains bound (not shown). In fact, as 
shown in Fig. 9 (see also Fig. 10), a second electronic 
state becomes localized about the hole above U ext ~ 
1.25 Ry*, corresponding to n e 0.2 a* B ~ , thus forming 
a so-called trion |26|. Here we have plotted the (Kohn- 
Sham) level energies for the localized hole case, measured 
relative to the band edge far from the hole, so that an en- 
ergy above zero indicates an asymptotically free particle. 



a = 0.35 (triangles) 0.45 (hexagons) and 0.55 a* B (boxes). We find that for even higher densities, > 8 x 10 6 
Line gives corresponding n e (see text for exact definition). 
Rise of energy below n e = caused by shift of electronic 
charge to gate for delocalized hole case. 



cm 



a third electron can become bound to the hole; a state 
which we refer to naturally as a "quadron." 



conduction band density n e , which is here defined as the 
total number of conduction band electrons in the wire for 
the localized hole case, minus unity (the bound electron) 
divided by the wire length L y . When the gate voltage be- 
comes negative, all electrons become depleted from the 
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FIG. 9. Discrete Kohn-Sham energy levels for electron 
(solid) and hole (dashed) vs. gate voltage, a = 0.35 a* B . To 
lowest order exciton binding energy is Eb ~ s\ + e\, hence 
binding weakens with increased Uext (and n e ). Hole localizes 
two electrons above U ex t ~ 1 -Ry*. 

wire and, below U ex t — — 1.0 Ry*, the conduction elec- 
trons for the delocalized hole case shift from the wire to 
the gate. This explains the jump in Eb below depletion. 
The minimum of Eb occurs at the lowest density before 




y (O 



FIG. 10. Occupied electronic wave functions (moduli 
squared), states 1 through 19, a = 0.35 a* B , n e — 0.6 a g ~ x . 
Two states are localized near hole (at y = 0). Orthogonal- 
ity of remaining states to bound states suppresses screening 
charge density near hole. 

The most intriguing feature of the results in Fig. 9 
is that, while the hole becomes more weakly bound, and 
consequently more spatially extended, the electrons, par- 
ticularly the lowest state, become more strongly bound 
with increasing U ex t (and consequently increasing n e ). A 
similar result for electron binding to an ionized impurity 
in a 2DEG silicon inversion layer was found by Vintcr 
p3[ . This is rather counter-intuitive, since one expects 
screening, by free electrons, of the interaction between 
hole and bound electron to weaken the attraction and 
separate the particles. Partially this result is understood 
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as relating to the direct product nature of the exciton 
state when expressed in DFT, as opposed to the single 
composite state of the two-body or screened two-body 
problem. Thus the binding energy is not a single eigen- 
value, or even the sum of two eigenvalues, but rather it 
must be understood as a difference between two inter- 
acting ground states. Nonetheless the decreased energy 
and increased localization of the lowest Kohn-Sham level 
seems puzzling. In figure 10 we show, for a relatively 
high density, n e = 0.6 a^ _1 , the set of all occupied eigen- 
functions (moduli squared). At this n e the lowest two 
electron states are localized near the hole at y = 0. Note, 
however, that due to the restraint of orthogonality, the 
densities of all other states in the vicinity of the hole are 
suppressed. Thus while the total conduction band 




-4 -2 2 4 

y (O 

FIG. 11. Total electron density (solid), total density minus 
bound state (dashed), effective electrostatic potential V e ff(y) 
(dotted), density of hole (dot-dash) and density of gate charge 
(double-dot-dash) versus y, same parameters as Fig. 10. 

density, exhibited in Fig. 11, peaks at y = 0, this den- 
sity minus that of the lowest state diminishes at the hole. 
Combined with the uniformly increasing background pos- 
itive charge on the gate, this comprises a build-up of pos- 
itive charge near the hole as U ex t and n e are increased. 
Therefore an increased density of electrons produces an 
effectively anti-screening effect which initially drives the 
electron closer to the hole and then admits a second and 
even a third bound electron. The hole eigenstate, mean- 
while, is growing increasingly de-localized and the faster 
rise of its energy accounts for the increase of the overall 
energy (i.e. a decrease of the negative binding energy). 

B. Combined results 

Finally, in combining the results from Section II and 
Section III, we simply add the excitonic binding energy 



to the bandgap calculated in LDA for the translation- 
ally invariant wire. The results for this are shown in 
Fig. 6. Clearly the exchange-correlation and kinematic 
effects causing the BGR produce a red shift which is 
partially negated by the shift of the exciton binding en- 
ergy due to screening. This is consistent with the no- 
table insensitivity of the photolumincscencc line position 
in optical experiments on quantum wires to either pho- 
toexcited electron-hole plasma or to gate-generated 
1DEG |l2|. Interestingly, it is also consistent with recent 
calculations of the dynamically screened Bethe-Salpeter 
equation by Das Sarma and Wang [Q for wires with a 
two-component plasma, although the interpretation dif- 
fers somewhat from that discussed here. While the ten- 
dency toward cancellation of BGR and excitonic screen- 
ing is concluded in both studies, the most striking differ- 
ence between the results of Ref . @ and those here concern 
the unbinding of the exciton which, in an electron-hole 
plasma, is called the Mott transition. The authors of 
Ref. [0 estimate a merging of the exciton energy with 
the electron- hole band edge at around n~3x 10 5 cm^ 1 , 
where n is here of course the density of electrons and 
holes. Even in the simple one-electron static screening 
approximation they find a vanishing of the exciton bind- 
ing energy with density. By contrast, we find that up to 
the highest density considered (3 x 10 6 cm -1 ) the exciton 
remains bound. Assuming that both conclusions are cor- 
rect, the implication is that it is the interaction between 
excitons which leads to their unbinding or, equivalently, 
their merging with the continuum. 

In any event, the continued existence, despite screen- 
ing, of a bound state of the electron with the hole in one 
dimension is not surprising in that in ID any small attrac- 
tive potential should bind an electron. One could expect 
that, at some density the screening of electrons near the 
Fermi surface would effectively cancel the total potential 
produced by the hole, thereby releasing the bound elec- 
tron. However for this to occur the Fermi wavelength 
would need to be much smaller than the exciton radius, 
which is ~ 1 a* B . Our calculations have not, as yet, shown 
a density where such a transition of the screening cloud 
from states at the bottom of the Fermi sea to states at the 
top of the Fermi sea occurs, but should such a transition 
exist it would be of considerable interest theoretically. 

IV. CONCLUSION 

In conclusion, we have presented results of density 
functional calculations for the electronic structure of 
modulation doped and gated T-shaped quantum wires for 
the case where the wire can be assumed to be translation- 
ally invariant. We have shown that the phenomenon of 
bandgap renormalization can be qualitatively and quan- 
titatively understood within DFT. We have further em- 
ployed a simplified model of a cylindrical quantum wire 
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to examine the strength of the bound exciton state in 
the presence of a one dimensional electron gas. We have 
normalized out the band structure problem here by defin- 
ing the exciton binding energy as the difference between 
the total energy of the wire with a localized hole and 
that with the hole charge (and consequently all elec- 
trons) spread uniformly along the wire. We find that 
the variation of the exciton binding energy with density 
tends to cancel the bandgap renormalization, in agree- 
ment with recent experiments. Finally, we have noted 
that, in this one component plasma case, to the highest 
densities which we have considered, no analogue of the 
Mott transition, i.e. no unbinding of the electron and 
hole occurs. Rather, the orthogonality of the free elec- 
tron states with those of the bound electron(s) leads to 
an anti-screening behavior such that as n e increases a 
second (trion) and even a third (quadron) bound state 
forms at the hole. 

In the future we hope to investigate the effect of dimen- 
sionality on the exciton physics by extending the calcu- 
lation to include multiple ID subbands. 
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